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S/P  TRAVEL  TIME  RATIO  IN  CENTRAL  ASIA 


1.  INTRODUCTION 


This  report  describes  work  accomplished  from  October  1,  1973,  through  31  March 
1974,  on  the  investigation  of  the  use  of  close-in  P  and  S  arrivals  to  deter¬ 
mine  focal  depth  of  seismic  events  in  Central  Asia.  Emphasis  of  the  study  to 
date  has  been  on  empirical  cal’bration  of  the  S/P  travel  time  ratio  for  earth¬ 
quakes  in  or  near  the  Hindu-Kush  intermediate  depth  source  observed  at  stations 
withm  13  degrees  of  the  source.  The  report  includes  an  evaluation  of  the  im¬ 
portance  of  depth  and  distance  dependence  of  the  S/P  travel  time  ratio. 

The  S/P  travel  time  ratio  for  shallow  Central  Asia  earthquakes  observed  at 
close-in  stations  is  currently  being  investigated. 

In  section  2  the  basic  background  required  to  understand  the  theory  behind  the 
method  used  is  briefly  discussed.  The  use  of  P  and  S  velocities  on  depth  of 
focus  and  origin  time  determinations  is  emphasized  together  with  the  calibra¬ 
tion  of  close-in  stations  and  the  use  of  ’’  ndati  graphs  (Wadati,  1933,  Kisslinger 
and  Engdahl ,  1973).  In  sections  3  and  4  the  data  used  in  the  study  and  the 
methods  employed  are  described.  In  section  5  the  results  obtained  are  discussed. 


2.  BACKGROUND 


The  modified  Geiger  method  (Macelwane  and  Sohon,  1932)  commonly  used  to  locate 
earthquakes  yields  more  accurate  focal  depth  if  origin  time  is  known  indepen¬ 
dently  of  the  P  arrival  times  used  to  mak®  the  location  (Evernden,  1969, 
Swanson,  1969,  1971,  Kisslir.grr  and  Engdahl,  1973).  Origin  time  (Ts)  can  be 
computed  from  S  and  P  arrival  times  at  a  single  station  if  the  ratio  of  S  to 
P  travel  times  is  known  for  the  path  from  hypocenter  to  station: 


=  P  -  (S 

-P)/r 

where  r 

=-  (vy  -1 

Ts 

=  S  travel  time 

TP 

=  P  travel  time 

S 

=  S  arrival  time 

P 

=  P  arrival  time 

TS 

=  origin  time  estimated  by  the  S/P  technique 
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Depth  fhs)  is  then  estimated  by  a  routine  location  procedure  with  origin  time 
restrained  to  Ts.  The  depth  accuracy  to  expect  when  locating  with  the  modi¬ 
fied  Geiger  technique  with  origin  time  restrained  may  be  clarified  by  the 
theory  of  hypocenter  control  factors  (Yeith  and  van  Leer,  1070) .  The  depth 
control  factor  for  a  given  station  is  defined  as  the  change  in  depth  estimate 
resulting  from  a  1  second  residual  at  the  station.  For  the  ideal  case  of  a 
symmetrical  net,  so  there  is  no  interaction  between  epicenter  and  depth  esti¬ 
mates,  the  depth  control  factor  for  station  i  is  given  by 


O/ah^ 

(depth  control  factor).  =  - y 

i  Efat/dhlj2 


The  change  in  depth  (dh)  is  given  by 


dh 


L(3t/3h)iUi 

i _ _ 

F.Ot/3h)i2 

i 


where  uj  =  travel  time  residual 
at  station  i_ 

If  origin  time  is  restrained  to  an  erro  eons  value  (T'  =  T  +  6T) ,  then 


dh' 


L(3t/3h)j  (u^-6T) 
i 


E  (nt/3h).2 

i  ^ 


E  (3t/3h). 

dh'  =  dh  +  — - - 

E(3t 

(6h)  resulting  from  an  error  in  origin  time  (5T) 
T(3t/3h)i 

J -  (-6T) 

z(;>t/3h).2 

i  1 
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Therefore,  the  error  in  depth 
is  given  by 

6h  = 


As  a  rule  of  thumb,  for  a  good  net  of  stations  (Swanson,  1969V. 


-  AT 

Ah  =  - 

(St  /  di)- 
A»  n 


where  A  =  mean  distance  of  locating  net. 


The  accuracy  with  which  r  (see  equation  1)  must  be  known  for  a  specified 
accuracy  in  Ts  may  be  estimated  easily,  assuming  no  arrival  time  errors  in 
S  or  P: 


Ar 


(S-P) 


AT 


where 

AT 

=  origin  time  error 

Ar 

3  error  in  r 

For  example,  assuming  a  value  of  0.78  for  r  and  an  S-P  time  of  25  seconds,  an 
error  in  r  of  .024,  or  one  part  in  33,  will  result  in  a  1  second  origin  time 
error,  while  for  an  S-P  time  of  100  seconds,  Ar  is  .006,  or  one  part  in  130, 
for  an  origin  time  error  of  1  second.  The  required  accuracy  in  r  for  a  speci¬ 
fied  origin  time  accuracy  increases  rapidly  with  increasing  distance. 

In  practice,  a  significant  portion  of  the  error  in  determining  the  S/P  travel 
time  ratio  results  from  errors  in  the  S  arrival  time,  either  because  of  diffi¬ 
culty  in  timing  the  start  of  the  phase  or  because  of  phase  misidentification. 

If  S  and  P  arrival  times  from  several  close-in  stations  are  available,  Wadati 
graph  analysis  (Kisslinger  and  Engdahl,  1973,  Wadati,  1933)  may  be  used  to 
determine  the  ratio  of  P  to  S  velocity  at  the  ray  bottom  and  the  origin  time, 
if  the  form  of  the  S-P  vs  P-0  relation  and  the  normal  P-0  intercept  value  are 
known.  In  the  most  general  case,  the  Wadati  graph  is  nonlinear  and  the  P-0 
intercept  is  not  zero.  Only  if  P  and  S  travel  times  are  linear  functions  of 
distance  is  S-P  a  linear  function  of  P-0;  however,  the  P-0  intercept  is,  in 
general,  still  not  zero. 

In  order  to  derive  the  relationship  let 

S-T0  =  ag  +  A/vs 

P-To  =  ap  ♦  A/Vp 

Then  S-P  =  [as  -  (p+l)  ap]  *  p  (P-T0) 

where  p  =  vp/vs  -  1 
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or 


S-P  =  p[P  -  (Tq  +  t)1  (2) 

whore  x  =  -  I  [as  -(p  +  1)  ap] 


Kaila  (1969")  and  Kaila,  Krishna,  and  N'arain  (1969)  estimated  P  and  S  velocities 
as  a  function  of  focal  depth  for  earthquakes  in  the  Hindu-Kush  intermediate 
depth  source,  using  time-distance  plots  for  distances  near  the  inflection  point 
of  the  time  distance  curves.  Because  the  ray  parameter  changes  only  slowly 
with  distance  over  this  range  of  distance,  the  ray  parameter  may  he  estimated 
from  the  slope  of  the  time-distance  curve  in  this  region.  Then  the  velocity 
at  the  hypocentral  depth  can  be  determined,  since  the  ray  corresponding  to  the 
inflection  point  on  the  travel  time  curve  leaves  the  hypocenter  horizontally 
(the  hypocenter  is  at  the  ray  bottom  for  this  path)  and  the  depth  is  presumed 
known . 

In  the  present  study,  IVadati  graph  analysis  has  been  carried  out  for  earth¬ 
quakes  in  the  Hindu-Kush  intermediate  depth  source,  using  limited  distance 
ranges  to  more  closely  approximate  linearity  of  the  IVadati  graph,  primarily  to 
investigate  depth  dependence  of  the  P-0  intercept.  Change  of  P-0  intercept 
with  depth  is  important  because  a  deviation  of  the  intercept  from  its  expected 
value,  resulting  from  the  difference  in  actual  and  initially  assumed  depths, 
looks  like  origin  time  error  when  using  the  S/P  technique  to  locate  events  of 
unknown  (or  uncertain)  origin  time. 


3.  EXPERIMENTAL  DATA 


P  and  S  arrival  times  from  43  intermediate  depth  Hindu-Kush  calibration  earth¬ 
quakes  recorded  at  6  IVhSSN  stations  and  19  other  close-in  stations  were  used 
in  the  study.  The  calibration  earthquakes  occurred  during  1964  through  1971. 
All  had  body  wave  magnitude  of  4.4  or  greater  and  had  previously  been  relocated 
using  preliminary  source-region/station  time  (SRST)  corrections  (Veith,  1973), 
and  with  focal  depth  restrained  to  the  depth  determined  from  pP-P  intervals. 
Herrin  (1968)  P  travel  times  were  used  in  making  the  locations.  These  locations 
are  considered  to  be  the  most  reliable  locations  available  for  the  calibration 
events.  Locations  of  the  calibration  earthquakes,  both  shallow  and  deep,  are 
shown  in  figure  1.  Table  1  gives  the  locations  and  origin  times  of  the  cali¬ 
bration  earthquakes.  Table  2  lists  station  codes,  locations,  and  elevations 
of  the  stations  whose  data  were  used  in  the  study.  In  this  table,  WWSSN  sta¬ 
tions  are  indicated  with  an  asterisk. 

P  and  S  times  of  the  events  recorded  at  the  six  WIVSSN  stations  were  reviewed 
from  film  reproductions  of  their  seismograms.  Arrival  times  were  changed  or 
added  when  appropriate.  P  and  S  arrival  times  at  the  other  close-in  stations 
were  taken  from  the  tri-monthly  seismological  bulletin  of  the  Moscow  Academy 
of  Science  and  the  bulletin  of  the  International  Seismological  Centre  (ISC). 

In  cases  where  times  in  Moscow  and  ISC  bulletins  disagreed,  times  were  taken 
from  the  Moscow  bulletin. 
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Figure  1.  Calibration  earthquakes 
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Table  1.  Calibration  earthquakes 


Date 
Yr.  Mo. 

Da. 

Origin  Time 

Hr:  Min:  Sec: 

North 

Latitude 

degrees 

East 

Longitude 

degrees 

hpp 

km 

nb 

71 

10 

14 

21:55:55.2 

36.51 

71.26 

86 

5.0 

67 

12 

17 

00:25:19.2 

36.52 

71.41 

97 

5.6 

70 

09 

05 

19:26:27.3 

37. 13 

71.42 

100 

5.0 

67 

08 

12 

22:54:39.4 

37.11 

71.33 

101 

5.1 

71 

06 

26 

22:23:31.4 

36. 38 

71.41 

121 

5.2 

70 

11 

13 

17:30:09.6 

36.98 

71.50 

124 

5.1 

65 

04 

10 

21:21:29.4 

37.44 

71.80 

124 

4.9 

68 

09 

27 

10:37:58.9 

37.88 

72.26 

124 

5.3 

67 

11 

07 

19:57:28.6 

37. 12 

71.77 

128 

5.0 

69 

11 

24 

17:23:23.0 

37.19 

71.66 

129 

5.8 

70 

01 

08 

21: 19: 13.7 

36.34 

70.89 

140 

4.9 

64 

12 

24 

01:08:39.3 

36.31 

70.94 

145 

5.2 

64 

03 

23 

13:40:31.1 

38.25 

73.61 

150 

5.2 

64 

03 

16 

03:28:16.2 

37.96 

72.83 

152 

5.0 

71 

03 

18 

19: 12:46.4 

38.22 

73.64 

153 

4.9 

68 

05 

08 

22:45:10.9 

37.17 

71.91 

161 

4.9 

71 

10 

15 

16:22:15.3 

36.99 

71.89 

165 

4.9 

70 

12 

26 

07:58:41.2 

36.53 

70.91 

165 

4.4 

71 

06 

22 

06:30:03.3 

36.29 

69.87 

167 

5.0 

64 

06 

06 

08:05:59.8 

37.13 

71.92 

171 

4.8 

70 

06 

11 

17:40:53.9 

36.47 

70.99 

195 

5.1 

70 

11 

29 

08:53:25. 1 

36.58 

71.36 

200 

4.8 

69 

09 

12 

05:08:04.3 

36.42 

70.91 

201 

4.9 

66 

08 

16 

02:16:22.9 

36.42 

70.79 

208 

5.5 

71 

08 

04 

00:24:39.5 

36.40 

70.69 

210 

5.6 

65 

10 

06 

15:35:07.2 

36.48 

70.11 

213 

5.0 

71 

08 

04 

01:59:06.3 

36.44 

70.73 

214 

4.9 

70 

07 

21 

01:18:08.0 

36.46 

70.36 

215 

5.1 

66 

03 

31 

23:38:04.3 

36.40 

70.71 

216 

5.4 

70 

09 

14 

16:11:18.6 

36.47 

70.14 

220 

4.8 

65 

03 

14 

15:53:09.2 

36 . 33 

70.71 

220 

6.4 

64 

11 

16 

04:47:29.9 

36.44 

70.54 

222 

5.  1 

67 

05 

08 

18:48:08.0 

36.47 

70.13 

223 

4.8 

66 

06 

06 

07:46:18.6 

36.39 

71.11 

225 

6.2 

70 

01 

26 

16:38:34.1 

36.46 

70.50 

228 

4.6 

66 

06 

04 

05:11:58.6 

36.35 

70.69 

228 

5.4 

69 

09 

04 

02:57:22.1 

36.51 

70 . 85 

231 

4.7 

69 

10 

28 

18:45:13.5 

36.31 

70.86 

234 

5.0 

65 

07 

24 

17:57:44.8 

36.44 

71.15 

238 

4.7 

65 

11 

16 

01:03:58.9 

36.38 

71.09 

252 

5.1 

65 

09 

14 

18:57:33.1 

36.39 

70.10 

261 

4.6 

70 

09 

04 

13: 12:03.7 

36.64 

70.17 

280 

4.7 

67 

01 

25 

01:50:22.2 

36.63 

71.57 

289 

5.7 
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Table  2.  Station  locations 


North  Hast 

Latitude  Longitude  Elevation 


Station 

Code 

deg 

min 

sec 

deg 

min 

sec 

meters 

Talgar,  Kazakh  SSR 

AAB 

43 

16 

00.0 

77 

23 

00.0 

850 

Andizhan,  Uzbek  SSR 

ANR 

40 

45 

18.0 

72 

21 

36.0 

494 

Ashkhabad,  Turkmen  SSR 

ASH 

37 

57 

00.0 

58 

21 

00.0 

220 

Dehra  Dun,  India 

DD1 

30 

19 

00.0 

78 

02 

00.0 

682 

Dushambe,  Tadzhik  SSR 

DSH 

38 

33 

30.0 

68 

46 

30.0 

817 

Fergana,  Uzbek  SSR 

FRG 

40 

23 

00.0 

71 

47 

00.0 

600 

Frunze,  Kirgiz  SSR 

FRU 

42 

50 

00.0 

74 

37 

00.0 

655 

Garm,  Tadzhik  SSR 

GAR 

39 

00 

00.0 

70 

19 

00.0 

1300 

Karachi,  Pakistan 

KAP. 

24 

56 

00.0 

67 

08 

36.0 

34 

Kizyl-Awat,  Turkmen  SSR 

KAT 

39 

12 

00.0 

56 

16 

00.0 

90 

Kabul,  Afghanistan 

KBL* 

34 

32 

27.0 

69 

02 

35.4 

19  20 

Khorog,  Tadzhik  SSR 

KUO 

37 

29 

00.0 

71 

32 

00.0 

1850 

kulyab,  Tadzhik  SSR 

KUL 

57 

54 

00.0 

69 

45 

00.0 

605 

Lahore,  Pakistan 

LAH* 

31 

33 

00.0 

74 

20 

00.0 

210 

Mangla,  Pakistan 

MNL 

53 

08 

50.0 

73 

45 

00.0 

436 

Meshed,  Iran 

MSH* 

30 

18 

40.0 

59 

55 

16.0 

987 

Murgab ,  Tadzhik  SSR 

MUR 

38 

22 

00.0 

73 

56 

oo.  n 

3800 

New  Delhi,  India 

ND1  * 

28 

41 

00.0 

77 

13 

00.0 

207 

Nilore,  Pakistan 

ML* 

35 

59 

00.0 

73 

15 

06.0 

536 

Narvn,  Kirgiz  SSR 

NRN 

41 

26 

00.0 

76 

00 

00.0 

2849 

Przhevalsk,  Kirgiz  SSR 

pr: 

42 

29 

00.0 

78 

24 

00.0 

1599 

Quetta,  Pakistan 

QUH* 

30 

11 

18.0 

66 

57 

00.0 

1721 

Samarkand,  Uzbek  SSR 

SAM 

39 

40 

24.0 

66 

59 

24.0 

704 

Tashkent,  Tadzhik  SSR 

TAS 

41 

19 

30.0 

69 

17 

42.0 

470 

Warsak,  Pakistan 

WRS 

34 

09 

00.0 

71 

25 

00.0 

345 

*WWSS.N  Station 


4.  METHODS 


1. 1  CALIBRATION  OF  CLOSE-IN  STATIONS 

lor  each  calibration  earthquake,  values  of  the  S/P  travel  time  ratio  r  [see 
equation  1"  were  computed  from 

r  =  (S-rmr-Tpp)  (3) 

where  T  p  =  origin  time  estimate  from  relo¬ 
cation  restrained  to  pP  depth 
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for  each  station  that  recorded  both  P  and  S  from  the  event.  Mean  values  of  r 
were  computed  using  the  station  mean  r  values  by 


6T 


lij  =  (p_TpP)ij  '  (s_p)ij/ri 


(4) 


where  r^  =  mean  r  for  station  ^ 

6Tiii  =  origin  time  deviation  for  station 
i_  and  event  j_. 

Station  and  event  origin  time  deviations  (S^  and  Ejj,  respectively)  were 
computed  using  the  station  mean  r  values  from 


n; 


su  ' 

1  3  =  1 


(5) 


where  nj  =  number  of  events  for  which 
station  i  had  data 


4,  ■  ^  £  "h, 

3  i=l 


(6) 


where  nj  =  number  of  stations  with  data 
for  event  ^ 

In  addition,  station  and  event  origin  time  deviations  (S2i  and  I^j* 
respectively)  were  computed  using  the  same  value  of  r  (the  overall  mean 
value)  for  all  stations: 


S2i  =  ~  L  1 

ni  3  =  1 


(7) 


E2  j 


II: 


h  £  [C’-Vij  -  1 

J  i=l 


(8) 


where  r  =  overall  mean  r 
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Results  are  given  in  section  5.3. 


4.2  WAD AT I  GRAPH  ANALYSTS 

To  invest' gate  the  sensitivity  of  the  S/P  travel  time  ratio  to  changes  in 
focal  depth,  Wadati  graph  anal>sis  was  performed  on  the  dat a  for  each  event, 
using  only  arrivals  within  a  -estricted  distance  range.  Based  on  the  results 
reported  by  Kaila,  and  others  (1969)  on  S  and  P  travel  times  from  Hindu-Kush 
earthquakes,  a  distance  range  of  4.5  to  13.0  degrees  was  selected,  to  assure 
that  the  assumption  of  a  linear  relation  between  S-P  and  P-0  is  justified. 
Reliable  estimates  of  Wadati  graph  parameters  could  not  be  obtained  for  in¬ 
dividual  events  because  too  few  stations  within  the  restricted  distance  range 
recorded  both  P  and  S  for  any  single  event.  Therefore,  data  from  earthquakes 
of  similar  depth  were  combined  and  least  squares  estimates  of  the  Wadati  graph 
slope  (p)  and  P-0  intercept  (t)  were  determined  for  each  of  10  depth  ranges. 
This  procedure  is  valid  if  the  pP  origin  times  are  correct  and  if  the  Wadati 
graph  parameters  do  not  change  rapidly  with  depth. 

The  Wadati  graph  for  intermediate  depth  earthquakes  should  not  be  linear  at 
distances  less  than  about  4.5  degrees.  However,  because  almost  half  the  data 
fall  in  this  distance  range,  it  is  of  interest  to  determine  if  S  and  P  arrival 
times  at  stations  in  this  distance  range  can  be  used  to  compute  accurate 
origin  times,  using  the  slope  and  P-0  intercept  determined  from  Wadati  graph 
analysis.  Wadati  graph  analysis  should  also  aid  in  investigating  sensitivity 
of  the  S/P  travel  time  ratio  to  changes  in  focal  depth  and  epicentral  distance 
over  this  distance  range.  Therefore,  linear  Wadati  graph  parameters  were  es¬ 
timated  for  each  of  the  10  depth  ranges,  using  only  arrivals  at  distances  less 
than  4.5  degrees. 

Using  the  weighted  mean  values  of  Wadati  graph  slope  (p)  and  P-0  intercept  (x), 
origin  time  deviations  from  the  pP  values  were  computed  by 


6T3ij  =  (P-TpP}ij  '  Tk  '  (S_P)ij/Pk 

where  k  =  distance  range  index 


(9) 


Mean  station  (S3i)  and  event  (E3O  origin  time  deviations  were  computed  by 
taking  the  appropriate  mean  of  the  deviations  defined  by  equation  (9).  The 
event  origin  time  deviations  were  computed  for  each  distance  range  separately 
and  for  the  two  distance  ranges  combined.  Then  residual  event  origin  time 
deviations  were  computed  by 


(10) 


where  nj 

Results  are  given  in  section  5.3. 


number  of  stations  with 
data  for  event  j. 


-10- 


TR  74-9 


5.  RESULTS  A\'D  DISCUSSION 


5.1  DEPTH  DEPENDENCE  OF  P/S  VELOCITY  RATIO 

Table  3  lists  the  Wadati  graph  parameters  estimated  from  the  data  grouped  by 
depth  range  for  the  two  distance  ranges,  4.5  to  13.0  degrees  and  less  than 
4.5  degrees.  Figure  2  shows  the  Wadati  graph  slopes  and  P-Tpp  intercepts, 
with  plus  or  minus  one  standard  error  limits,  as  a  function  of  pP  depth  for 
the  distance  range  4.5-13.0  degrees,  with  the  mean  slope  and  intercept  indi¬ 
cated.  Also  shown  on  this  figure  are  values  of  Wadati  graph  slope  (pi  and 
P-0  intercept  (t)  computed  from  the  P  and  S  velocity  values  reported  by  Kaila, 
et  al  (1969)  for  the  liindu-Kush  region  and  smooth  curves  of  p  and  i  (see  equa¬ 
tion  2)  as  functions  of  hpP  computed  from  Kaila' s  linear  velocity-depth  func¬ 
tions.  Most  of  the  difference  between  Kaila' s  intercept  values  and  the  inter¬ 
cept  values  determined  in  the  present  study  can  be  accounted  for  by  the  dif¬ 
ferent  P  travel  times  used,  Jeffreys -Bui len  (195ft!  (JR)  in  Kaila's  study  and 
Herrin  (1968)  in  the  present  study.  Locations  made  with  Herrin  times  yield 
origin  times  about  2.5  seconds  later,  on  the  average,  than  locations  made  with 
JB  times. 

The  minima  in  the  P  to  S  velocity  ratio  (p+1)  and  P-0  intercept  in  the  140  to 
170  km  depth  ranges  appear  to  be  real.  However,  in  view  of  the  error  limits 
on  slope  and  intercept  estimates  for  the  individual  depth  ranges,  constant 
slope  and  intercept  values  over  the  full  depth  range  constitute  a  plausible 
interpretation  of  the  data.  Deviations  of  the  intercept  values  from  the  over¬ 
all  mean  value  (+  1.6  second)  range  from  +2.1  to  -1.9  seconds.  Tf  interpreted 
as  origin  time  errors,  and  assuming  a  mean  depth  derivative  of  about  -0.1  sec¬ 
ond/km  for  events  in  this  depth  range,  corresponding  depth  errors  would  be 
about  ±20  km. 

Wadati  graph  parameters  as  a  function  of  pP  depth  for  distances  less  than  4.5 
degrees,  with  the  mean  slope  and  intercept  indicated,  are  shown  in  figure  3. 

For  this  distance  range,  rays  do  not  leave  the  focus  horizontally.  Therefore, 
the  Wadati  graph  slope  can  not  be  interpreted  in  terms  of  the  ratio  of  P  to  S 
velocity  at  the  hypocentral  depth.  However,  the  technique  provides  a  means 
of  using  data  from  stations  in  this  distance  range  to  compute  origin  times 
from  P  and  S  arrival  times,  recognizing  that  the  assumption  of  a  linear  Wadati 
graph  is  only  approximate  in  this  case. 

5.2  DISTANCE  AND  DEPTH  DEPENDENCE  OF  THE  S/P  TRAVEL  TIME  RATTO 

Figures  4  through  8  show  observed  values  of  r  (  =  (S-P) /(P-TpPl)  as  a  function 
of  P-Tpp  for  each  of  the  10  depth  ranges  used  in  the  Wadati  graph  analysis. 

The  smootn  curves  shown  or  each  figure  represent  the  expected  value  of  r,  com¬ 
puted  from  the  mean  Wadati  graph  parameters  for  the  two  distance  ranges. 

The  large  increase  in  observed  r  values  at  very  close  distances  is  accounted 
for  by  the  strong  negative  P-0  intercept  for  the  Wadati  graph  for  this  distance 
range  (figure  3).  For  distances  greater  than  4.5  degrees,  distance  dependency 
of  the  S/P  travel  time  ratio  is  n*.  strong,  the  slight  increase  in  r  with  in¬ 
creasing  distance  being  accounted  for  by  the  positive  Wadati  graph  intercept 
for  this  distance  range  (figure  2).  The  strongest  evidence  of  depth  dependency 
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rs  for  distances  less 


travel-time  ratio  for  depth  ranges  36-101  km  and  121-129  km 


Sure  7-  S/P  travel-time  ratio  for  depth  ranges  220- 


of  r  occurs  for  the  161  to  171  km  depth  range,  for  which  the  curve  computed 
from  the  mean  Wadati  graph  parameters  gives  r  values  mostly  larger  than  ob¬ 
served  values  and  the  252  to  261  km  depth  range,  for  which  the  r  values  com¬ 
puted  from  the  mean  Wadati  graph  parameters  are  smaller  than  observed  value 
(see  also  figures  2  and  3).  In  general  as  expected  from  results  of  the  Wadati 
graph  analysis  (figures  2  and  51,  the  curves  computed  from  the  mean  Wadati 
graph  parameters  represent  the  observed  S/P  travel  time  ratios  better  for  le 
4,5  to  13  degree  distance  range  than  for  the  less  than  4.5  degree  range. 


5.3  S/P  ORIGIN  TIME  ESTIMATES 

Principal  results  of  the  study  to  date  are  summarized  in  table  4,  giving  sta¬ 
tion  statistics  and  table  5,  giving  statistics  on  estimates  of  event  origin 
time  deviations,  using  the  different  techniques  described  in  section  4. 

Table  4  lists  the  mean  r  and  standard  deviation  of  r  for  each  station,  the. 
mean  and  standard  deviation  of  origin  time  deviations  from  the  p  value  using 
the  station  mean  r  for  all  events  (Sul,  using  the  same  value  ((.780)  of  for 
all  stations  (S2i) ,  and  using  the  mean  Wadati  graph  parameters  (S31)  .  Table  5 
summarizes  event  ori jin  time  deviations  from  the  pP  value  computed  with  the 
various  techniques:  using  the  station  mean  r  s  (E^,  equation  6)  both  using 
all  stations  and  using  only  stations  whose  data  were  reviewed  (with  the  re¬ 
quirement  that  at  least  two  stations  have  da*a  for  each  event!,  using  the 
same  value  of  r  (0.780)  for  all  stations  (E2j,  equation  8),  and  using  the  mean 
Wadati  graph  slope  and  intercept  for  the  appropriate  distance  range,  for  the 
two  distance  ranges  both  separately  and  combined.  In  addition,  the  residual 
event  origin  time  deviations  (E4j,  equation  10),  after  removing  the  station 
mean  P-0  intercepts  (S5i)  are  included  in  table  5. 

Figure  0  shows  the  event  origin  time  deviations,  computed  by  the  different 
techniques ,  as  a  function  of  pP  depth.  No  clear  trend  with  depth  is  apparent 
for  orig.n  time  deviations  computed  by  any  of  the  four  techniques,  indicating 
that,  for  identification  purposes,  depth  dependency  of  the  S/P  travel  time 
ratio  or  the  Wadati  graph  parameters  is  not  serious  for  this  region.  Locations 
of  the  two  earthquakes  with  the  largest  residual  origin  time  deviations 
(Ei 4: +2 . 5  seconds  at  165  km  and  -2.2  seconds  at  261  km)  were  based  on  pP  depths 
less  well  determined  than  those  of  the  other  calibration  earthquakes;  there¬ 
fore,  errors  in  the  calibration  depths  may  account  for  these  large  deviations. 


As  is  evident  from  the  results  shown  in  table  5,  there  was  little  difference 
in  origin  time  deviations  from  the  pP  values  computed  by  the  different  tech- 
niques  based  on  S/P  travel  time  ratios.  Using  each  station  s  mean  r  (E^)  did 
as  well  as  using  the  estimates  of  slope  (0  =  Vn/vs  -  1)  and  intercept  (t)  from 
the  Wadati  graph  analysis  (E3O  with  the  two  distance  ranges  combined,  both 
techniques  giving  a  standard  deviation  of  origin  time  deviation  )f  0.9  second, 
lor  events  in  the  intermediate  depth  range,  the  corresponding  deoth  error 
(deviation  from  pP  depth)  would  have  a  standard  deviation  of  8-9  km. 


Removing  the  mean  station  effect  (Sji)  resulted  in  a  further  slight  reduction 
in  standard  deviation  of  origin  time  deviation  from  the  pP  value  (E4j  m  table 
5).  Because  of  the  large  values  of  the  mean  station  effects  for  some  stations 
rg,!  in  table  4)  it  is  important  to  correct  for  these  factors  when  using  on  v 
a  few  stations  to  compute  origin  times  by  means  of  the  Wadati  graph  parameters. 


-20- 


TR  74-9 


T  ab  1  e  4 . 

Station  S/P  travel  time 

ratios 

and 

origin  time  deviations 

No .  of 
events 

Mean 

dist. 

(deg) 

r 

mean 

St  dv 

secs 

Sli 

mean  St  dv 

secs 

$2i 

mean  st  dv 

secs 

S3i 

mean  st 

32 

1.18 

.807 

.047 

0.2 

1.5 

-0.8 

1.5 

0.2 

1. 

22 

1.75 

.813 

.037 

0.2 

1.4 

-1.3 

1.5 

-0.6 

1 . 

9 

2.32 

.772 

.013 

0.0 

0.7 

0.4 

0.7 

0.4 

0. 

28 

2.50 

.780 

.051 

0.2 

2.6 

0.1 

2.6 

0.2 

7 

2.73 

.778 

.029 

0.1 

1.5 

0.2 

1.5 

-0.2 

1 

36 

2.74 

.778 

.038 

0.0 

2.0 

0.2 

2.0 

0.0 

2 

11 

2.93 

.758 

.045 

0.2 

2.8 

1.5 

2.8 

0.9 

2 

9 

3.56 

.782 

.020 

0.2 

1.3 

0.0 

1.3 

-1.2 

1 

8 

3.58 

.778 

.010 

0.1 

0.7 

0.2 

0.7 

-1.3 

0 

29 

4.05 

.773 

.025 

0.0 

1.9 

0.5 

1.9 

-0.8 

2 

13 

4.14 

.758 

.032 

0.0 

2.7 

1.8 

2.7 

-0.2 

2 

12 

4.38 

.737 

.036 

-0.1 

3.1 

3.5 

2.9 

2.0 

2 

33 

4.84 

.771 

.020 

0.1 

1.7 

0.9 

1.8 

0.0 

1 

26 

5.89 

.760 

.014 

0.0 

1.6 

2.1 

1.5 

1.5 

1 

5 

6.04 

.792 

.015 

-0.1 

1.6 

-1.4 

1.7 

-1.9 

1 

27 

6.71 

.776 

.015 

-0.1 

1.8 

0.4 

1.8 

0.0 

1 

30 

7.35 

.784 

.011 

0.0 

1.4 

-0.5 

1.4 

-0.8 

1 

15 

8.36 

.792 

.014 

0.0 

2.0 

-1.8 

2.0 

-1.9 

2 

13 

8.49 

.785 

.031 

0.0 

4.6 

-0.8 

4.7 

-0.9 

4 

15 

8.58 

.781 

.013 

0.0 

1.8 

-0.2 

1.9 

-0.3 

1 

20 

9.45 

.783 

.013 

0.0 

2.1 

-0.5 

2.1 

-0.5 

2 

WWSSN  station  -  data  reviewed 

- 
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Table  4  (Continued) 


Mean 

se 

cs 

secs 

secs 

St  a- 

No.  of 

dist . 

r 

S 

li 

S2i 

^5i 

tion 

events 

(deg) 

mean 

St  dv 

mean 

St  dv 

mean  St  dv 

mean  St  d\ 

NDI* 

33 

9.35 

.767 

.011 

0.1 

1.8 

2.1  1.8 

2.1  1.8 

ASH 

12 

10.33 

.797 

.015 

0.1 

2.7 

-2.9  2.7 

-2.7  2.7 

KAT 

13 

12.02 

.  ”89 

-0.1 

3.4 

-1.9  3.4 

-1.5  3.4 

KAR 

10 

.771 

.017 

0.1 

3.7 

2.0  3.7 

2.4  3.6 

*WWSSN  station  -  data  reviewed 


Table  5.  Event  origin  time  deviations 


Event 

origin  time  P-0 

Station 

Mean 

St  dv 

deviation  intercept 

Slope 

correction 

secs 

secs 

m 

1—* 
t— 1. 

o 

n 

0 

All  stas 
reviewed  stas 

-0.01 

0.86 

and  2  or  more 

1.07 

E2j  0 

r 

0 

All  stas 
reviewed  stas 

+0.25 

and  2  or  more 

1.17 

E3j  ~k 

ck 

0 

both  A  ranges 

0.90 

A  <  4.5° 

-0.09 

1.08 

4.5°  <  A  <  13.0° 

+0.02 

1 . 09 

E4j 

s3i 

both  A  ranges 

0.81 

A  <  4.5° 

1.21 

4.5°  1  A  £  13.0° 

0.91 

i  :  station  index 

j  :  event  index 

k  :  distance  range 

index 
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Figure  9 


hpp  (km) 


Origin  time  deviations  from  TpP 
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For  identification  purposes,  using  the  same  value  of  r  for  all  stations  Cl^2 j ^ 
is  entirely  adequate  for  Hindu-Kush  intermediate  depth  earthquakes,  the  stan¬ 
dard  deviation  of  origin  time  deviations  for  this  case  being  1.0  second,  with 
a  corresponding  standard  deviation  of  depth  error  of  about  10-32  km. 


6.  CONCLUSION'S  AND  RECOMMENDATIONS 


The  S/P  travel  time  ratio  technique  of  computing  focal  depth  will  identify 
llindu-Kush  intermediate  depth  earthquakes  if  reliable  S  and  P  arrival  times 
at  stations  within  13  degrees  of  the  epicenter  are  available.  Variation  of 
the  S/P  travel  time  ratio  with  epicentral  distance  is  not  significant  for  the 
distance  range  4.5  to  13  degrees.  For  distances  less  than  4.5  degrees,  the 
S/P  travel  time  ratio  increases  rapidly  with  decreasing  distance.  For  the 
llindu-Kush  intermediate  depth  source,  depth  dependence  of  the  S/P  travel  time 
ratio  is  weak  and  can  be  neglected  for  identification  purposes. 

We  recommend  that  the  investigation  of  the  S/P  travel  time  ratio  in  Central 
Asia  be  continued.  Close-in  stations  should  be  calibrated  for  the  S/P  travel 
time  ratio  for  shallow  earthquakes  in  this  region.  Sensitivity  of  the  ob¬ 
served  travel  time  ratio  to  changes  in  focal  depth  should  be  investigated. 

The  depth  accuracy  obtainable  for  shallow  Central  Asia  earthquakes  should  be 
determined. 
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